4e-condensation in a fully frustrated Josephson junction diamond chain 
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Fully frustrated one-dimensional diamond Josephson chains have been shown [B. Dougot and J. 
Vidal, Phys. Rev. Lett. 88, 227005 (2002)] to posses a remarkable property: The superfluid phase 
occurs through the condensation of pairs of Cooper pairs. By means of Monte Carlo simulations 
we analyze quantitatively the Insulator to 4e-Superfiuid transition. We determine the location of 
the critical point and discuss the behaviour of the phase-phase correlators. For comparison we also 
present the case of a diamond chain at zero and 1 /3 frustration where the standard 2e-condensation 
is observed. 



PACS numbers: 

Josephson arrays in the quantum regime have been 
studied extensively 0, both experimentally and theo- 
retically, as model systems where to investigate a variety 
of quantum phase transitions. The application of a mag- 
netic field creates frustration and leads to a number of 
interesting physical effects 0, 

Very recently, renewed interest in frustrated Joseph- 
son networks has been stimulated by the work by Vidal 
et al. 13] on the existence of localization in fully frus- 
trated tight binding models with symmetry. Local- 
ization in this case is due the destructive interference for 
paths circumventing every plaquctte. These clusters over 
which localization takes place were named Aharonov- 
Bohm (AB) cages. Experiments in superconducting net- 
works have been performed and the existence of the AB 
cages has been confirmed through critical current mea- 
surements both in wire Q and junction |^ networks. 
Starting from the original paper by Vidal et al. several 
aspects of the AB cages both for classical 0, H IS ^ 
and quantum [lll | superconducting networks have been 
highlighted. 

The basic mechanism leading to the AB cages is 
also present in the (simpler) quasi-one-dimensional lat- 
tice shown in FigC] At fully frustration, it has been 
shown that superconducting coherence is established 
throughout the system by means of 4e-condcnsation. The 
global superconducting state is due to the condensation 
of pairs of Cooper pairs. Predictions on the critical 
current of the diamond chain of Fig^ amenable of ex- 
perimental confirmation have been put forward by Pro- 
topopov and Feigelman 0,0. Unusual transport prop- 
erties of these systems have been also predicted in semi- 
conducting samples |l5l |. 

In this work we present the results of our Monte Carlo 
simulations on the Josephson junction network with the 
geometry depicted in Fig^ Our aim was to perform a 
detailed quantitative analysis of the phase diagram pre- 
dicted in Ref.^3l- order to have a fairly complete 
description of the effect of frustration in this case, we 




FIG. 1: The diamond chain Josephson network analyzed in 
the present paper. The crosses represent the Josephson junc- 
tions connecting two neighboring superconducting islands. In 
the chain there are two types of inequivalent sites with two 
(grey) and four (black) neighbours. By an appropriate choice 
of the gauge, the magnetic phase factors Ai^j can be chosen 
to be zero on the three links indicated by a continuous lines 
and / in the fourth one indicated by a dotted line. 

considered the stiffness and phase correlators for three 
values of the frustration parameter; i.e. / = 0, / = 1/3 
and / = 1/2. 

The Hamiltonian for a Josephson junction network is 
7^-So5I"?-^^IIcos(^, . (1) 

The first term in the Hamiltonian is due to the charg- 
ing energy. Here for simplicity we consider the case in 
which the Coulomb interaction is on-site; see Ref.fisj for 
the more realistic case of long range charging interac- 
tion. The second term is the Josephson contribution. 
The phase of the superconducting order parameter in the 
i-th island is denoted by ipi , Eq is the charging energy and 
Ej is the Josephson coupling energy. The number and 
phase Lpi operators are canonically conjugate on each site 
[ui , 6**^^ ] = 5ij e**^' . The gauge-invariant definition of the 
phase in presence of an external vector potential A and 
flux-per-plaquette <I> ($o = hc/2e is the flux quantum) 
contains the term At , = ^ f.'' A • dl . All the obscrvables 

*U $0 J'l 

are function of the frustration parameter defined as 
f = i- I A • dl A,, 
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where the hne integral is performed over the elementary 
plaquette. Due to the periodicity of the model it is suf- 
ficient to consider values of the frustration < / < 1/2. 

The Monte Carlo simulations have been performed on 
an effective classical d + 1-dimensional XY- model 0, 0| 
(here d = 1) whose action is given by 

S = -K ^ cos {ifii^k - <fiiM') 
i,{kk') 

-K ^ cos {ipi^k - ^j,k - Ai^j) . (2) 

The effective dimensionless coupling is defined as K = 
Ej I Eo . The index k labels the extra (imaginary time) 
direction which takes into account the quantum fluctua- 
tions (the vector potential does not depend on the imagi- 
nary time). The first term corresponds to charging while 
the second is due to the Joscphson coupling. The simu- 
lations where performed on L x L lattice with periodic 
boundary conditions (the largest lattice was 72 x 72). In 
Eq.(|21) the couplings along the time and space directions 
have been made equal by a proper choice of the Trotter 
time slice 0, ^3 ■ This choice, with no consequences on 
the study of the zero temperature phase transition, makes 
the analysis of the Monte Carlo data considerably sim- 
pler. The expectation values of the different observables 
(stiffness and correlation functions) have been obtained 
averaging up to 10'' Monte Carlo configurations by using 
a standard Metropolis algorithm. Typically the first half 
of configurations, in each run, were used for thcrmaliza- 
tion. 

The stiffness, related to the critical current, is used to 
signal the presence of the transition. It is defined through 
the increase of the free energy f due to a phase twist S 
imposed along the space direction ■ 



The critical point is expected to be of the Berezinskii- 
Kosterlitz-Thouless universality class 0, l^^l . Its loca- 
tion can be determined using the following ansatz for the 
size dependence of T{Kc) 



rUKc) = 1 



1 



21n(i/;o) 



(3) 



where lo is the only fit parameter. In the presence of 
frustration, the universality class of the transition may 
be different from that of the unfrustratcd case. In the 
case of the two-dimensional fully frustrated XY-model 
this issue has been investigated in great detail (see Ref. 
I21I and refs. therein). Up to date, there is no unanimous 
consensus on the nature of the transition. However, in 
this work we suppose that the transition belongs to the 
BKT universality class, as suggested by Ref.^3|j ^.nd de- 
termine the critical value by means of Eq. Q . 
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FIG. 2: In the upper panel, the stiffness for the case of / = 
frustration is plotted against the coupling. Different symbols 
correspond to different sizes of the chain: circles to L = 36, 
squares to L = 48, diamonds to L = 56, triangles up to 
L — 64, and triangles down to L — 72. The dashed line with 
slope 2/tt gives a rough estimate of the transition point. A 
better estimate is obtained by means of the finite size scaling 
shown in the lower panel and explained in the text. Thick 
black line has exactly slope 1 and is plotted as a reference 
guide. The value of lo at the transition is 2.9 



We first analyze the / = case and extract the value 
of the critical coupling from the stiffness. This extrapo- 
lation has been done by performing a linear fit in loga- 
rithmic scale [tt KTl{K) — 2] ^ ~ InL — \nlo and 
searching for the coupling value such that a(K) = 1. This 
coupling value is then identified with the critical point 
Kc- The proposed ansatz fits very well the data and the 
estimated value of the critical coupling is K^^ = 1.28 
which corresponds to {Ej/Eq)c ^ 0.61. Data are re- 
ported in Fig|21 

The results of the stiffness at / = 1/3 and 1/2 are 
shown in Fig|31 and Fig0| respectively. As compared to 
the unfrustrated case, the critical value of the Joseph- 
son coupling required to establish superfiuid coherence 
is slightly larger for / = 1/3 and further increases for 
the fully frustrated case / = 1/2. The ansatz of Eq. Q 
seems to provide an accurate estimate of the transition 
point for / = and / = 1/3. In the fully frustrated case, 
however, the value of = 25 indicates that we probably 
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FIG. 3: The same plots of FigHfor the case of / = 1/3. The 
critical point is K^^ = 1.045 (with lo ~ 0.6). 

need larger chains in order to really enter the critical re- 
gion. Another indication of this fact emerges in the upper 
panel of Fig0]where the line of slope 2/7r crosses the data 
already when the stiffness is decreasing to zero. In order 
to put bounds to the critical point in the fully frustrated 
case we plot in Fig|31the stiffness clS 0j function of the 
system size. From the raw data it is possible to bound 
the transition point in the range 0.55 < < 0.57. 
All these results are summarized in table below: 



/ 





1/3 


1/2 


{Ej/Eo)c 


1.28 ±0.01 
0.61 ±0.01 


1.045 ±0.005 
0.91 ±0.01 


0.56 ±0.01 
3.2 ±0.1 



The ratio of the obtained critical couplings for unfrus- 
trated and fully frustrated system is i^c, 1/2/^0,0 = 2.28± 
0.06 and not 4 as expected from the reduction by a factor 
1/2 of the effective charge of the topological excitation 
that unbind at the critical point. This may be due to 
the fact that the screening of the vortices is different in 
the unfrustrated and fully frustrated case therefore lead- 
ing to a further correction in the ratio between the two 
critical points. 

The differences in the fully frustrated case manifest 
dramatically in the way condensation is achieved. As pre- 
dicted by DouQot and Vidal j the destructive interfer- 
ence built in the diamond structure prevents Cooper pair 
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FIG. 4: The same plots of FigsEEfor / = 1/2. Compared 
to the cases of / = and / — 1/3, the superfluid region is 
considerably shrunk. The value of lo at the transition is ~ 25. 

to have (quasi-)long range order. The superfluid phase is 
then established via the delocalization of pairs of Cooper 
pairs. This is at the origin of the 4e— condensation. In or- 
der to check this point, the knowledge of the phase-phase 
correlators is required. Quasi-long range behaviour in a 
two-point correlation function of the type 

g2n{\i - j\) = (cosn((^, - (pj)) (4) 

signals the existence of condensation of 2n charged ob- 
jects. In Figl^lwe discuss their properties. In the upper 
panels, we consider the phase-phase correlator 32 for two 
different couplings deep in the superfluid and Mott in- 
sulating phases respectively. What is evident from the 
figure is that, despite the fact that the system is phase 
coherent, phase correlations decay very fast almost in- 
dependently from the value oi K. As explained in p^ . 
this behaviour should be ascribed to the existence of the 
Aharonov-Bohm cages. Even if hopping of single Cooper 
pairs is forbidden because of quantum interference, cor- 
related hopping of two pairs does not suffer the same 
destructive interference. In the lower panels of the same 
figure, the space dependence of the correlator 174 is plot- 
ted for the same coupling as upwards. The different be- 
haviour between the Mott and the superfluid phase is 
now evident. The correlator decays exponentially only 
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FIG. 5: The stiffness is plotted as a function of tlie system 
size for different values of K in the critical region. This plot 
highlights the existence of a transition, though do not allow 
to extract the transition point. For K"^ < 0.55 data seem 
to scale to a finite value in the thermodynamic limit, whereas 
over 0.57 it seems clear that they go to zero. Different symbols 
correspond to values of K~^: circles (0.45), squares (0.49), 
diamond (0.53), triangles up (0.555), triangles left (0.56), tri- 
angles down (0.565), triangles right (0.575). 




FIG. 6: The phase correlators (?2 and are shown as a 
function of the distance for the fully frustrated case / = 1/2. 
Data are plotted for a chain with L = 48. On the left side 
circles correspond to = 0.1 deep in the ordered phase 

and squares to — 0.5 on the border of it. On the right 

side, squares are = 0.6 and circles = 1.2, deep in 

the Mott insulator phase. Differently from g2, the correlator 
Qi shows quasi- long range order. 



for K^^ = 1.2 > Kc (right side): in the other panel, 
differently from 172, the decay is power-law like. For com- 
parison we report also simulations of the phase correla- 
tors for the case f = 1/3. In this case the "standard" 
condensation of Cooper pair is observed as witnessed by 
the behaviour of 52 shown in FigEI 

We acknowledge useful discussions with Michclc Gov- 




FIG. 7: Phase correlators at frustration / — 1/3. Up: the 
phase correlator g2{\i — j\) is shown as a function of the dis- 
tance between the sites both in the ordered phase (left panel) , 
— 0.3 (circles) and 1.0 (squares) and in the Mott insula- 
tor phase (on the right) = 1.1 (squares) and 1.4 (circles). 
Differently from the fully frustrated case, here the phase cor- 
relator of Cooper pairs changes its behaviour at the critical 
point. Down: the phase correlator g4,{\i — j\) is shown for the 
same coupling values as upwards. 
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